7th  Int.  Symp.  “Nanostructures:  Physics  and  Technology’ 
St  Petersburg,  Russia,  June  14-18,  1999 
© 1999  Ioffe  Institute 


NT.  07 


Self-organized  growth  of  composition-modulated  alloys 

V.  A.  Shchukin  and  A.  N.  Starodubtsev 

Ioffe  Physico-Technical  Institute,  St  Petersburg,  Russia 


Epitaxial  growth  of  semiconductor  alloys  in  open  systems  can  be  unstable  against  fluctua- 
tions of  composition  [1-4],  For  such  kinetic  instability  it  is  typical  that  the  bulk  diffusion 
is  negligibly  slow  compared  with  the  surface  diffusion.  Fluctuations  of  composition  in 
the  entire  thickness  of  the  epitaxial  film  are  “frozen”,  and  these  fluctuations  create  a strain 
field  on  the  surface  which  affects  the  surface  migration  of  atoms.  It  was  shown  in  [4], 
that  due  to  the  interplay  of  anisotropic  elasticity  and  anisotropic  surface  diffusion,  the 
wave  vector  of  the  most  unstable  mode  of  composition  fluctuations  can  have  any  direction. 
However,  an  important  question  was  still  unsolved:  How  does  this  interplay  depend  on  key 
experimental  parameters,  namely  on  growth  temperature  and  growth  velocity?  Moreover, 
the  linear  stability  analysis  of  alloy  growth  [1-4]  does  not  give  an  answer  on  the  final 
composition-modulated  structure  which  is  being  formed  during  the  growth. 

In  the  present  paper  we,  first,  perform  a linear  stability  analysis  and  obtain  the  wave 
vector  of  the  most  unstable  mode  as  a function  of  growth  temperature  and  growth  velocity. 
Second,  we  solve  a non-linear  problem  in  the  weak  segregation  regime  and  obtain  a final 
composition-modulated  structure. 

We  study  the  growth  of  a binary  alloy  A|_,;Bt.  via  molecular  beam  epitaxy  at  growth 
temperature  above  the  roughening  transition  temperature.  Then  the  surface  contains  a high 
concentration  of  steps  and  kinks.  Crystal  growth  can  be  described  by  the  propagation  of  the 
surface  in  the  normal  direction.  The  growth  of  an  alloy  implies  coupled  fluctuations  of  alloy 
composition  z)  — c(i*||,  z)  — c and  surface  profile  £(r||)  = h( ry)  — h,  where  h(r\\) 
is  the  local  thickness  of  the  film,  and  h is  the  average  thickness.  Let  the  alloy  A]_CBC  be 
a mixture  of  two  cubic  crystals,  the  homogeneous  alloy  be  lattice-matched  to  the  substrate 
and  the  growth  occur  on  the  (00 1 (-substrate  of  a cubic  material.  The  growth  of  an  alloy  is 
described  by  two  coupled  kinetic  equations.  Let  the  reference  frame  travel  together  with 
the  surface  with  the  average  growth  velocity  v.  Then  kinetic  equations  are  as  follows: 
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Here  Mfj  and  M?.  are  2D  surface  atomic  mobility  tensors  related  to  the  substitutional 
diffusion  and  to  the  change  of  surface  profile,  respectively.  Unlike  [3],  Eqs.(l)  account  the 
anisotropy  of  surface  mobility  tensors,  their  principle  axes  being  [110]  and  1 1 10],  which  is 
typical  for  III-V  and  II-V1  semiconductors;  a is  the  lattice  parameter.  The  free  energy  F 
equals  Fchem  + Fgrad  + Esurf  + Eei,  where  Fchem  is  the  chemical  free  energy,  Fgrad  is  the 
gradient  energy,  Esurf  is  the  surface  energy  of  a corrugated  surface,  and  Eej  is  the  elastic 
energy  depending  on  the  surface  profile. 

First,  we  carry  out  the  linear  stability  analysis  by  expanding  F up  to  <fi2  and  £2.  We 
seek  Fourier  components  of  fluctuations  0(k|| , t)  and  l (kp , t)  proportional  to  explcu  (ky  )t ]. 
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Fig.  1.  The  phase  diagram  of  kinetic  instability  in  variables  “temperature-growth  velocity”, 
(a)  The  boundary  of  the  instability  region.  For  v = vo  the  alloy  is  unstable  at  T2  < T < T\ . (b)  I: 
stable  alloys,  II:  unstable  alloys,  most  unstable  modes  are  shown  in  Fig.  1(c)  (bottom).  Ill:  unstable 
alloys,  most  unstable  modes  are  shown  in  Fig.  1 (c)  (top),  (c)  Wave  vectors  of  most  unstable  modes. 


Linearized  equations  for  0(ky , t)  and  f (kg,  t)  are  decoupled,  similar  to  [2,3 J-  For  fluc- 
tuations of  composition,  0(k|| , t)  ~ exp[cu(k|)f],  we  obtain  the  following  characteristic 
equation  for  o>: 


co  = —M(0)k2 
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0 is  the  angle  between  ky  and  [100];  Rs(0)  refer  to  the  anisotropic  elastic  interaction  [1  ]; 
as(0)  are  attenuation  coefficients  of  static  Rayleigh  waves;  M (6 ) is  the  effective  mobility 
in  the  given  direction,  M{0)  = cos2  (0  — tt  / 4)  + sin2 (0  — n/4),  where  Mf  and  M? 

are  mobilities  along  [1 10]  and  [1 10],  Tc  is  the  critical  temperature  of  kinetic  instability  in 
the  slow  deposition  limit  (u  -*■  0)  [1];  r = (93  fci,em/dc2dT),  where  fchem  is  the  chemical 
free  energy  density;  k refers  to  the  gradient  energy.  Given  v and  T,  Eq.(2)  determines 
co(k,  0).  Retu  > 0 means  the  kinetic  instability. 

To  obtain  the  stability  phase  diagram  of  the  alloy  growth  (Fig.  1 (b))  we  set  c = 0.5,  use 
for  fchem  the  regular  solution  approximation,  fchem(c,  T)  — Qcfl  — c)  + T [ ( 1 — c)  ln(l  — 
c ) + c In  c],  and  take  into  account  the  Arrhenius-type  temperature  dependence  of  mobilities 
M\'2{T)  — m\°U*P(-EU2/T).  The  growth  of  an  alloy  includes  the  competition  between 
the  surface  migration  and  the  burial  by  the  incoming  flux.  Let  temperature  T be  fixed.  Then, 
if  the  growth  velocity  v is  larger  than  a certain  critical  value,  the  fast  burial  does  not  allow  the 
surface  migration  to  create  composition  fluctuations.  Now,  let  v be  fixed  and  T decrease. 
When  T < T\  (v),  the  driving  force  to  kinetic  instability  appears.  When  T < IMv)  the 
slow  surface  diffusion  does  not  allow  composition  fluctuations  to  occur  (Fig.  1(a)). 

Anisotropic  surface  diffusion  promotes  instability  with  k|  along  the  direction  of  fast 
diffusion  [110].  Anisotropic  elasticity  promotes  instability  with  kj  along  the  elastically 
soft  direction  [100]  or  [010].  The  result  of  such  an  interplay  is  given  in  Fig.  1(b).  In  region 
II  there  are  2 most  unstable  modes  with  k;  ||  [ 1 10].  In  region  III  there  are  4 most  unstable 
modes  with  ky  deviated  by  ± 80  from  [1 10]  (Fig.  1(c)). 

Below  we  seek  steady  state  solutions  in  the  weak  segregation  regime,  i.e.  at  T and  v 
close  to  the  onset  of  the  instability.  For  simplicity  we  consider  isotropic  surface  diffusion. 
Near  the  onset  of  instability  we  have  ( VK)/\aMr(Tc  — T)\  — 1/4  — e2,  where  e <SC  1. 
The  fact  that  in  the  linear  regime,  kinetic  equations  for  $ and  for  < are  decoupled  and  the 


NT.07 


499 


surface  remains  flat,  has  two  important  issues.  First,  f ~ cp2.  Then,  to  solve  the  non-linear 
problem  we  expand  the  free  energy  up  to  fourth  order  terms,  i.e.  up  to  cp 4 and  f 2.  Second, 
the  relaxation  of  the  surface  is  a fast  process,  while  the  relaxation  of  the  composition 
profile  is  a slow  process.  By  using  the  adiabatic  elimination  of  fast  relaxing  variables  [5] 
we  obtain  kinetic  equations  which  contain  only  cp.  In  the  weak  segregation  regime  the 
solution  is  a perturbation  series  in  e.  The  main  contribution  is  given  by  most  unstable 
modes  with  ky  = ±(&o,  0)  and  ky  = ±(0,  ko),  where  ko  — (1  /\/2)[ r(Tc  — T)/k]]j2.  We 
seek  cp  in  the  form  cp(x,  y)  — €[(r^/r(Tc  — T)\ll2[\px  cos(kox)  + \py  cosf&oy)],  where 
r4  = ( 34  fchem / 9 c4 )c=c • Kinetic  equations  for  amplitudes  \px  and  \py  are 
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r = ( Mr(Tc  — T)/K)t  is  a dimensionless  time  parameter,  and  y is  the  surface  energy 
per  unit  area,  Rl  and  Rj  are  characteristic  elastic  energies.  The  steady  state  solution  of 

Eqs.(3)  depends  on  parameters  cil  and  ay.  If  ul  < My.  the  stable  steady  state  solution 

— 1/2  —1/2 

is  a ID  composition-modulated  structure,  \j/x  = aL  and  i/y  = 0 or  \f/y  — aL 
and  \j/x  — 0.  If  cil  > My\  the  stable  steady  state  solution  is  a 2D  modulated  structure, 
\j/x  — xfry  — ( ml  + cit  V1/2-  If  m/.  + ay  < 0 or  cil  < 0.  Eqs.(3)  do  not  have  stable  steady 
state  solutions.  This  means  that  the  final  structure  has  a finite  modulation  amplitude  already 
near  the  onset  of  instability. 

The  pattern  selection  is  governed  by  the  interplay  of  elastic  and  surface  energies.  Let 
the  surface  energy  be  very  high,  y — > oo.  Then  the  surface  remains  flat  even  in  the 
non-linear  regime;  cil  — L ay  — 3,  which  favors  a ID  structure.  There  is  a similarity 
with  alloy  decomposition  in  the  bulk  [6],  where  the  resulting  structure  is  a mixture  of  two 
phases  with  compositions  c\  and  a ■ A two-phase  system  is  consistent  with  ID  structure 
and  inconsistent  with  2D  structure.  For  a moderate  y,  the  effect  of  surface  corrugation 
becomes  important,  which  alters  cil  and  ay  and  favors  2D  structure.  This  is  similar  to  the 
case  of  lattice-mismatched  systems  where  a 2D  corrugation  gives  a more  efficient  elastic 
relaxation  than  a ID  one  [7],  For  a low  surface  energy,  m/,  < 0,  ay  < 0,  and  the  final 
structure  has  a finite  amplitude  near  the  onset  of  instability.  Since  the  interplay  of  elastic  and 
surface  energies  depends  on  Tc  — T (Eq.(4)).  this  yields  a phase  diagram  on  Fig.  2.  Similar 
interplay  between  ID  and  2D  structures  persists  in  case  of  anisotropic  surface  mobility. 

To  conclude,  we  have  shown  that  the  alloy  growth  in  open  systems  can  result  in  both  1 D 
and  2D  modulated  structure  and  the  orientation  of  the  structure  is  determined  by  growth 
temperature  and  growth  velocity.  Principle  results,  obtained  for  binary  alloys,  remain  valid 
for  ternary  and  quaternary  alloys.  This  explains  the  ID  structure  in  InAlAs  modulated  in 
[110]-direction  [8],  the  ID  structure  in  MgZnSSe  modulated  in  [100] -direction  [9],  and  the 
2D  structure  in  GalnAsP  modulated  in  [100]  and  [010]-directions  [10]. 
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Fig.  2.  The  phase  diagram  of  final  composition-modulated  structure  near  the  onset  of  instability. 
I:  ID  structure,  II:  2D  structure.  III:  no  solutions  in  the  weak  segregation  regime,  the  structure  has 
a finite  amplitude. 
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